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SUMMARY 


This report is intended to provide sufficient instructions for interfacing 
the Mangler-Smith, leading-edge vortex-rollup program with a vortex lattice 
(POTFAN) method and an advanced higher order, singularity linear analysis for 
computing the vortex effects for simple canard-wing combinations. 


INTRODUCTION 


Lifting surfaces generate vortices that can strongly Influence their own 
aerodynamic characteristics and those of nearby surfaces. Because many air- 
craft fighter configurations permit and even attempt to take advantage of such 
vortex-surface interactions, it is important to be able to predict the effects 
of these interactions. A technique has been developed that combines an approx- 
imate solution for the leading-edge separation and vortex rollup off highly 
swept lifting surfaces with two linear-panel methods to predict the vortex 
interactions. 

This manual is intended to demonstrate how to interface this leading-edge, 
vortex-rollup program with a vortex-lattice (POTFAN) code and an advanced 
higher order linear analysis referred to as the advanced panel code (APC) . A 
brief outline of the leading-edge separation, vortex-rollup analysis, and 
method of solution are also presented. 


MANGLER-SMITH METHOD 


The leading-edige separation and vortex-rollup analysis presented here is 
used to obtain a geometrical definition of the separated sheet and subsequent 
vortex rollup. By using this geometrical definition, we modeled the vortex 
sheet in the linear analyses. 

The theory (refs. 1-3) assumes slender-body conical flow in the cross- 
flow plane, giving a two-dimensional Laplacian as the governing equation. This 
means that the flow field at each cross section is independent of the flow at 
any other station and that variations in the stream direction x are much 
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smaller than those In the other directions, at least near the body. Therefore, 
the full potential <t> satisfies the equation 

(1 - + 4 . - 0 (1) 

XX yy zz 

where the first terra is negligible for highly swept wings. The general solu- 
tion of equation (1) is of the form 

it(x,y,z) - U (cos a) + U (sin a) + (})(x,y,z) (2) 

X z 

where U is the free-streain velocity and a is the angle of attack. The 
perturbation potential and its derivatives vanish at infinity. 

Combining equations (1) and (2) shows that tj> satisfies the two dimen- 
sional Laplaclan 

4“ 41 “0 (3) 

^yy -zz 

The conical-flow assumptions also mean that x must appear as a parameter in 
41 , which must therefore be expressible in the form 

♦ (x.y. 2 ) - . i) (4) 

The equation of the free vortex sheet in cylindrical polar coordinates 
(fig. 1 (a)) is assumed to be of the form 

S(x,r,e) Hr- sf(e) - 0 (5) 

where s is the wing semispan at some x station. So that 

s ■ X cot A 

where A is the leading-edge sweep angle. Two boundary conditions are applied 
to the sheet: 

1. The vortex sheet is a stream surface, so that 

(U cos a + 41 )S + 4> S + 4' S - 0 ( 6 ) 

^x X r r n n 

Assuming |4>^| << U, equation ( 6 ) in conjunction with equation (5), becomes 

4 , - ^ sin X (7) 

n X 

2. There is a zero pressure jump across the sheet, so that 

AC « 0 (8) 
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Using Bernoulli's equation 
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and therefore 




( 10 ) 


but 

^ ♦ 2) . ^(^ 2 + .^ 2) . 2(* ) (10a) 

y z on ^ m ° 

where the subscript m denotes the mean tangential velocity across the sheet. 
Again, from equation (4) 




± - Z * - * * 

X X ^y X 

1 A(j) - ~ A(y(*iy + i A4> - r cos XA 


(10b) 

(10c) 


By substituting equations (10a) and (10b) In equation (1) , the pressure condi- 
tion reduces to 


A(^ ■ A(^^|r cos X 




(U> 


Figure 1(b) shows the cross-flow plane Z ■ y + iz at some x-statlon. 

The free sheet ABC from the leading edge A terminates at C, where the spiral 
rollup begins. For purposes of estimating the leading-edge wake shape, all 
the voT.'tlclty within the spiral can be represented by an Isolated vortex of 
strength 2uyUaR/x, situated at the center D of a small circle of radius R, 
which joins smoothly the free sheet at G. Using equations (7) and (11), we 
can show that the equation of the spiral as It approaches D Is given by 

r' - # (12) 


where (r', 6') are polar coordinates relative to D. 

The calculations are carried out In a transformed plane Z* with the 
transformation (fig. 1(c)) 


2*2 « - s2 (13) 

in which the vortex sheet is assumed to be represented by a circular arc 
^*B*C* of radius Ci, and length 2r\6, and the rollup of the fed sheet is 
assumed to be eqivivalent to a small circle of radius r 2 and center D*. The 
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sheets join smoothly at C*. The transformation equation (13) enforces the 
requirements that the flow must be everywhere parallel to the wing surface. 

It is assumed that the vorticity distribution on A*B'*'C* is given by 
y(^) ' Yq ^r ^ “ ''^i 


The seven parameters (ri, r^, 0, u, y^» Y^* Y^) are determined from the fol- 
owing seven equations derived from the boundary conditions and mathematical 
compatibility requirements. 

1. Because the strengths at of the sheet and its image in the imag- 

inary axis must be equal and of opposite sign, then 

Y « 0 at A*(c - -0) 

Y^ + Yj. cos 0 + Yj! sin 0*0 (15) 

2. The free sheet A*B*C* must join the spiralling fed sheet smoothly 
at C*, so that the strength is continuous at C* (C ■ 0) 

2nUr2 /ai\^ 

Y + Y cos 0 - Y. sin G * { — ) (16) 

'o'r i x\a/ 

where the right side is the strength of the fed sheet at C*, aj being the 
distance A*D*. 

3. The singular point of the transformation equation (13) is A*; there- 
fore A* must be a stagnation point. Mathematically, this is equivalent to 
the condition 

2pr2a2 Yr Y^ 

180 “ ajs ~ ttU ® ® ttU ® 0) ■ 0 (17) 

where a 2 is the distance from D* to its image on the 0*Z* axis. 

4. There is a zero resultant force condition on the vortex system simi- 

lar to that in reference 4. This is equivalent to the condition that the non- 
singular part of the Induced complex velocity at D must be KU/s(2Z_ - Z ). 
This results in two equations: ^ 

5. the normal velocity condition 


applied at and 


-rU 


sin X 
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6. The zero pressure Jump condition 

^ - A(^^|r cos ^ ^ 

applied at 

The variables Yi» eliminated from equations (15)-(17) to give four 

nonlinear equations for r^, rzt 0* Yq with 

0<r2<ti; O<0< n/2 (18) 

These are solved through a least-squares minimization using the CONMIN optimi- 
zation package (ref. S). 


STREAMLINE GENERATOR MODULE 


The prediction of vortex rollup and vortex surface Interactions for com- 
plete aircraft configurations must be made numerically. However, numerical 
methods require an Initial estimate for the shape of the shed vortex sheet. 

In certain cases (Including in particular, the case of rollup from the leading 
edge) , the initial estimate must be reasonably accurate to guarantee conver- 
gence of the numerical method In a reasonable amount of time. Furthermore, a 
reasonable first guess might permit some useful results to be computed by 
using the trajectories of the vortices as input to POTFAN or APC. The Mangler- 
Smith theory provides such an initial estimate for the shape of the leading 
edge separated sheet. To compute the flow field, the sheet is discretized so 
that the separated wake is essentially modeled by a concentration of line vor- 
tices, which are described using vortex quadrilaterals. The vortex descrip- 
tion is thp'.i combined with a panel code in which the wakes are represented by 
doublet panels, which are either constant in strength (POTFAN) or have a qua- 
dratic variation (APC). To align the vortex trajectories with the stream- 
lines, the corner points along the leading-edge wake are obtained by integrat- 
ing the velocity field. 

The streamline generator module used is capable of producing trajectories 
of vortices shed from the leading, trailing, and side edges of lifting sur- 
faces. Input consists of POTFAN geometry files describing the lifting surfaces 
and data necessary for carrying out the least-squares minimization as described 
by equations (15)-(18). Output consists mainly of POTFAN wake geometry files 
for interfacing the corner points defining the wake trajectories with POTFAN 
and AFC. The streamline generator aligns the leading-edge shed wakes with the 
streamlines by integrating the flow field 

^ - V - (V • N)N (19) 

where V is the resultant velocity as computed by the Mangler-Smith theory 
and is the normal to the vortex-sheet equation (5) 
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N - vs 




( 20 ) 


where ^ Is a unit velocity along the scream direction and (r,,t\) are the 
two-dlml^slonal position and normal to Che trace of the sheet In the cross- 
flow plane. 

For the free sheet the normal Is along Che radius of the circular arc In 
the Z* plane. The corresponding components In Che Z plane are given by 
the transformation equation (13). For the spiral-fed sheet, the equation In 
the cross-flow plane Is given by equation (13) so that 

f(r'.e') i r’ - fr - 0 



( 21 ) 


Equation (19) was Integrated using an Improved fourth order Runge-Kutta scheme. 


POTFAN STRUCTURE AND ORGANIZATION 


The basic operational structure for the POTFAN system Is shown in 
figure 2 and consists of a batch of programs written In modular form. Each 
of the modules Is in an Independent main program by Itself. The modules com- 
municate with each other through files. The control statements required to 
handle the files are supplied automatically by Che program through the mecha- 
nism of command-format programming, which uses words or acronyms called 
"Commands" to control the action taken by the program. In addition, the pro- 
gram:. handles problems of variable size using dynamic memory allocation, thus 
enabling It to cope with large problems using a fixed amount of memory. 

The operations are Initiated through the POTGEM module that creates the 
geometry files for the lifting surfaces. POTWAK uses the Information stored 
on these files to create a wake file for the corner points defining the vortex 
trajectories, VVIM then generates the Influence matrices for the lifting sur- 
faces while INFMAN adds In the effect of the semi-independent wakes to create 
a flow field v. The boundary conditions on the lifting surfaces are calcu- 
lated in BCDN. PSOLVE calculates the singularity strengths by solving the 
system of linear equations of the form 


^ • N - -U • N 
(U + v)N^ ■ 0 


( 22 ) 
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Finally, the forces and monents are calculated In POTFOR. Further 
details on some of the POTFAM modules are given In references 6 and 7. To 
run the program requires the following Job control cards: 

POTGEM' create geometry file 

POTWAK* create wake file 

WIM> calculate Influence matrices 

INFMAM* add In semi-independent wakes 

BCDN> calculate boundary conditions 

PSOLVE* calculate singularity strengths 

POTFOR* calculate and print forces and moments 

7/8/9 

Each card must begin in column 1. Any Information typed in after the period 
Is considered a comment and may fill the entire card. The corresponding data 
Input decks are 

7/8/9 POTGEM DATA. 


7/8/9 POTWAK DATA. 


7/8/9 POTFOR DATA. 

In column 1 7/8/9 is typed, leaving column 2 blank. The organization of the 
data required for each of the POTFAN modules Is described In the following 
sections. 


POTGEM MODULE 


This module generates panel comer-point descriptions for the lifting 
surfaces. It involves the specification of suitable coordinate systems for 
characterizing the lifting surfaces together with automatic panel generation 
schemes. Basically, the Input to POTGEM consists of "commands," which arf: 
usually followed by data. The coimoands direct the logic flow and the data 
describes curves, panel distributions, singularity types, body rotations, ref- 
erence lengths, etc. The output consists of printed output and one or more 
geometry files. The geometry files contain the panel corner points, boundary 
condition points, unit normals, and sundry other data. The geometry files are 
used by the other POTFAN modules that require geometry or singularity data. 

Some typical input commands and their general ordering are described 
below: 

TITLE 

DSEGMENTS - defines panels on each network 


CARYySRIl “■ geometry definition 
SL,SU,VL,VU — network boundary definition 
SLBC,VLBC corner and control point distribution definition 
GRID ~ prints distributions defined In 5. 

RUSS rotate, shift and scale factors for network. 

PANL - 

DSFLAG type of shed vortex f quadrilateral, horseshoe, etc.) 

UVW — reads In unit wake vecttWs 
FINISH - 

ST^x^.E stores the current geometry on a disk file or tape. 

PRINT — prints geometry file 
STOP 

A number of specific examples are given In the following sections Illustrating 
the sample cases that have been analyzed. The specific Input requiring change 
to run different but similar configurations Is Identified. 

Table 1 Illustrates the Input deck for a flat delta wing of Al > 1.5 and 
A ■ 74°, with a root chord 3.4874 and semlspan 1.0 (fig. 3). To run a specific 
case, only the following >:ards need to be changed: 

Card 7 — NBPS ■ number of panels In the S-dlrectlon (chordwlse) 

NBPV ■ number of panels In the V-dlrectlon (spanwlse) 

Table 1 shows that, with ray paneling, there are 16 panels In the chordwlse 
and 8 panels In the spanwlse directions. For purposes of simplicity, a uni- 
form panel spacing Is used. More general lattice arrangements are possible 
and the Interested reader is referred to reference 6 for the corresponding 
Input layouts. 

Card 18 — VAR2.SV ■ root chord 

Card 20 — Along the leading edge VL referring to figure 3, VARISV speci- 
fies the limits within which x varies, while VAR2SV specifies 
the bounds on the y variable. 

Card 38 — FLT(l) ■ reference area for normalizing forces and moments. 

FLT(2,3,4) * reference lengths for normalizing moments about 
Oxf Oy, O 2 (fig. 3) axes respectively. 

FLT(5,6,7) ■ direction cosines of the semi-infinite straight 
wakes. If these wakes leave the wing plane at an angle a/ 2, 
these numbers become cos a/2, 0, and sin a/2. 

In the example, the wakes leave the wing parallel to Its plane. Therefore, 
the direction cosines are equal to 1, 0, and 0. 

Card 40 — Contains the geometry file Identifier that may be changed at 
the user's discretion. 

All Integers are Input under a 1615 format with floating point numbers under a 
8F10.1 format. In addition, all command formats must begin In column 1 (for 
example, cards 24-36, of table 1 and cards 37-57 of table 2). 
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Canacd-Ulng Comblnat^iOAi 


Tablas 2 and 3 contain tha corraspondlng input dacka for tha ainpla 
canard-wing combination of figura 4 and raportad in refaranca 8. Tha canard 
and wing ara aaaumad coplanar. Tha wing aamiapan ia 25.4 cm with a maan gao- 
matrical chord of 23.31 cm and araa 516.1 cm^. Tha following carda only naad 
ba changad in tabla 2 for tha canard gaomatry dafinition with atandard panala 
(flga. 4 and 5): 

Card 5 NBPS ■ numbar of panala in tha S dlractlon (apanwiaa) . 

NBPV ■ Numbar of panala in tha V diraction (chordwlaa) , 

In the example there ara 10 panala apanwiaa and 5 chordwlaa. Note that tha 

S, V directiona ara interchanged (figa. 3 and 5). 

Carda 19. 21 ~ tha variable VAR2SV daflnaa tha range of variation of tha 
y coordinate of tha leading edge. 

Carda 23, 24 - along tha leading edge VL» VARISV ■ range within which 
tha y-i;oordlnata varlaa, while VAR2SV apaclfiea the 
correat^oMding valuea for tha x variable. 

Card 26 — along tha trailing edge VU, VAR2SV <■ range within which 
X variaa 

Cards 34. 35 — RSHIFT(2) ■ y-eoordinata of tha old origin in the new 
system obtained by suppressing the fuselage. 

Cards 53, 56 — the diraction cosines of the semi-infinite wakiss from the 
trailing and 'side edges. 

In table 2, these leave at an angle of 20** to the canard plane. 

Card 59 ~ FLT(1,2,3,4) ■ as defined in table 1. 

FLT(IO) - X coordinate of point about which moments are taken. 

The default value is 0 in table 1. 

A similar setup prevails in table 3 for the wing, where the cards needing 
change are: 

Cards 5, 19, 21, 23, 24, 26 -• changes analogous to those in table 2 for 
the canard. 

Carda 39, 40 — corresponds to cards 34 and 35 of table 2 . 

Card 45 — this defines the point on the wing leading edge at which 

separation begins. This Is defined In terms of the cor- 
responding spanwlse panel on wing. 

Carda 58, 61 — corresponds to cards 53, 56, of table 2. 

Card 64 — corresponds to card 59 of table 2. 

In table 3, separation begins from the 11th panel. This la slightly outboard 
of the canard. If the panel number were 10^ then separation from the wing 
would begin at the point of Intersection of the canard tip line with the wing 
leading edge. 
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Card 32 — the AVLC command adjuata the corner and grid-point deacriptlona 
of the Ini'nraectlona of the wing leading edge with the corner- 
end control-point grid linea. The required nondlmenaional 
corner-point location XADJ ia calculated from the formula 


^corrected y coordinate of canard tii 
V wing aemiapan 


Thia number ahould lie between (-1, 1). The corner point nunr- 
est the XADJ will be made equal to XADJ. 


To calculate the corrected y coordinate of canard tip, POTGEM ia run on a 
stand-alone baais for the canard. A partial liating of the correaponding 
POTGEM output is shown in table 4 under the heading PANEL CORNER POINTS. The 
suffixes I, J denote the spanwise and chordwise variations respectively, 
with J ■ 1, corresponding to the leading edge of the canard. The 
y-coordlnate of the last panel point <I ■ 11) is the required point. Its 
y coordinate is 13.1121951, so that 


XADJ - 2 


’ 13.1121951 \ 
,25.40 - 3.81/ 


0.5 


0,214654479 


As an added check, the corresponding y coordinates up to I " 11 for the 
POTGEM wing program must be identical for the wing and canard. The spanwise 
panels for the canard and wing were generally selected according to the rule 


( NBPS)c ^ canard semispan 
(NBPS)w wing semispan 


POTWAK Module 

This module generates the corner points for the vortex trajectories by 
solving the equations of the Mangler-Smith theory. The leading-edge wake- 
corner points are generated by an Improved Runge-Kutta numerical integration 
of the flow field equations 19. In addition, there is a graphics capability 
for the display of the vortex trajectories. The program expects the following 
data to be entered by the namelist $DATA. 

NSEGS — number of segments for discretizing each wake line from the lead- 
ing edge. For highly swept wings at moderate angles of attack a 
value of 150 is generally found to be adequate. As the angle of 
attack decreases, NSEGS must be increased. 

ALPHA — angle of at' ack in degrees. Default is 15". 

UINF -- free-stream velocity. Default Is 1.0 in nondlmenslonalized units. 

ID ~ Identification number of the corresponding lifting surface geometry 
file. For the wing-alone case this number is given On card 40 of 
table 1. 
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PLOT ~ Loglciil variable whoae truth cauaas the waka and lifting surface 
geonetry to be plotted on the E&S picture system. Default is .FALSE. 

STORE — logical variable whose truth causes the wake and lifting surface 
geometry information to be stored on permanent file. Default is 
.TRUE. 


IDWAKE *- identification number of the wake file. IDWAKE is not required 
if STORE-. FALSE. 

DELTAX ~ X direction step used in calculating the leading edge vortex 
trajectories. Default of 0.1 in nondimensionallzed units. DELTAX 
must Increase in proportion to NSEGS. 

R1 -* estimate of radius of the circular arc in the Z* plane represent- 
ing the free sheet. Default in nondimensionalized units is 0.0765. 

R2 estimate of radius of the circle representing the fed sheet in the 
Z* plane. Default is 9.25E-3. 

THETA — estimate of half angle subtended by the finite circular arc in 
the Z*> plane. Default is 1.24 rad. 

GAMO — estimate of parameter associated with the vorticlty distribution 
on the free sheet. D/sfault in nondimensionalized units is -0.06384. 

'^TASHN ~ station code for the disposal of the plot file. The initial 

default value is 1, wtilch will cause the plot file to be sent to 

the FAE PDP-11. A value of 2 causes the file to be sent to the 

RJE terminal in Building 227 and a value of 3 causes it to be sent 

to the central-site printer. STASHM is not required if PLOT-. FALSE. 

IRKN — to reduce NSEGS to manageable levels, the integration scheme uses 
a number of intermediate points to generate a single Runge-Kutta 
point. This number iL represented by IRKN. For angles of attack 
that are sufl ciently small either NSEGS must be large or IRKN must 
be Increased. This also applies if the wing Is not sufficiently 
swept. Suggested values are NSEGS - 150 and IRKN - 10. 

VLB, VUB, IPRINT, NCON, NSIDE, INFO, UMAX, ICNDIR, NSCAL, SCAL, ITRM, 

CTL, CTLMIN, PHE, DELFUN, DABFUN, FDCH, FDCHM, THE, NFDG, LINOBJ, CT, 

CTMIN — CONMIN parameters defined in reference 5. The CONMIN parameters 
PHI and THETA are denoted here by PHE and THE. For the optimization 
technique to converge it was generally found that the finite differ- 
encing parameters FDCH and FDCHM for estimating the function deriva- 
tives must be small ~l.E-9. This is because of the presence of 
trigonometric functions in equations (15)-(17), which could give 
rise to multiple solutions for the least-squares minimization. 

These parameters need not be changed in normal use. 
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IFIAG* ILINEi lYEO, lYHI, IZLO, 12111, lARGMIN, lARUMAX - pAr«r.»&tev8 used 
in the subteutlnes lor the picture system es defined in reference 9. 
These parameters need not be changed. 

A typical listing la shown In table 5 Cor the delta wing whose configura- 
tion Is described in table 1, corresponding to an angle of attack a « 20", 

The delault values for the Mangler-Smith parameters Rl, R2, THETA, and GAMO 
are those appropriate to a 60" delta wing at 15" angle of attack, These 
default values were found to be adequate for all the cases run and gave very 
rapid convergence. For a wing with ray paneling, the number of wakos shed 
from the leading edge is given by 

NLINES - NRPS + 1 


Currently, NSEGSi*«NLlNES 10000, 

The pictui*e system displays the leading-edge vortex trajectories in the 
cross-flow plane, the (y,a) coordinates being scaled with respect to the 
stream coordinate x (fig. 6). The conlcul-<low assumptions imply that these 
trajectories must coincide for different x-stations. The straight segments 
represent the semi-infinite line vortices from the trailing edge, This is 
followed by orthogonal perspectives of the trajectories, figure 7, 

Table 6 shows the organ' latlon of the wake geometry output files, con- 
sisting of a number of records arranged according to the FOTFAN fortnats. 

Record 1 contains information relating to the flow-field characteristics, the 
number of wakes NLINES shed from the leading edge, the number of points 
MSECS describing each wake, the wake file identification, and other sundry 
POTGEM data. Record 2 contains the wake corner-point arrays XWAKE, YWAKE, 
and 2WAKE. Records 3-5 contain miscellaneous information that enables FOTWAK 
to converse smoothly with the other POTFAN modules as and when the need arises. 

To run the canard-wing case it is necessary to call POTGEM twice with 
different identifications. For example, the canard file could have an 
ID “ 8701, while the wing has an ID « 8702. This is then followed by two 
calls of POTWAR with corresponding lD*e in its input. The output files from 
POTCIEM and POTWAR can now be used to access the other POTFAN modules needed 
for the aerodynamic calculation. 


VVIM Module 

The VVIM module computes the influence matrices for the lifting surfaces, 
each element of the matrix representing the effect of a unit singularity such 
as horseshoe or quadrilateral vortices. VVIM is needed to form part of the 
set of simultaneous linear equations that represent the boundery conditions 
appropriate for given angles of attack. Input to VVIM coneiats of geometry 
files and user input. The user input consists of commands (which direct the 
logic flow) and data needed to effect the commands, The major commands in 
VVIM are the following: BREAD commands, which cause the geometry of sending 
(influencing) and/or receiving (influenced) components to be read inj IMAGE, 
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which causes the appropriate Iniage matrices and offset vectors to be calcu- 
lated; COMPRESSIBILITY » which causes geometries to be stretched to account for 
compressibility effects and which causes the velocity stretching data to be 
computed; SETWAKES, which can be used to change the wake shedding direction; 
and INFl, which causes the actual Influence matrix calculation and storing. 

Output from WIM consists of printed output and Influence matrix files. 
The Influence matrix files are used by PSOLVE, which Is the equation-solving 
module, and by POTFOR, which Is used to calculate forces, pressures, etc. 

Table 7 shows the WIM Input deck for the wing-alone case of tables 1 
and 3. No changes are needed In this deck to run different cases. Card 4 
shows that the POTGEM file Is accessed through ID » 30701. Card 6 instructs 
WIM to store the boundary condition Information on a disk file with 
ID * 31^01. NRWMAX Is the number of lifting surface networks. 

Table 8 shows the WIM Input deck for the canard-wing case of tables 2 
and 3. This table is divided Into four blocks; each block Is similar In 
structure to table 7. This division corresponds to a partitioning of the 
Influence coefficient matrix for the lifting surface Into four submatrices 
representing the canard-wing Interactions 



The caiiard-wlng case was run at a Mach number of 0.3 (card 15), so that the 
Prandtl-Glauert compressibility effects must be Included in the calculations, 
while table 7 corresponds to a zero Mach number, so that the COMP command Is 
not needed. The canard POTGEM file was accessed through ID * 8701, while the 
ID for the wing geometry file Is 8702. The resulting partitions Aq^, 

and are stored on disk file with respective ID's ■ 8801, 8802, 

8804, and 8803. The only card which needs to be changed Is card 15 for the 
Mach number. To calculate the singularity strengths T, it Is necessary to 
augment the lifting surface. Influence matrix A to include the effects of 
the wake networks. This Is done In the INFMAN module. The output from WIM 
consists primarily of the normal components of the induced velocities due to 
the lifting surfaces. In INFMAN the effects of the shed wakes are added to 
the normal velocity components to give v.N term In equation (22). 


INFMAN Module 

Wakes may be treated either as completely independent networks or as 
parts of the lifting surfaces to which they are attached. Treating wakes as 
Independent networks Is disadvantageous because unnecessary unknowns are Intro- 
duced and because pressures on panels adjacent to wake attachment lines will 
be incorrect. Treating wakes as Integral parts of the lifting surfaces to 
which they are attached is disadvantageous because this limits wake shapes to 
be in the form of singly ruled surfaces, which is unsatisfactory in the case 
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o£ leading- and side-edge rollup, or both, because It results In Inaccuracies 
during calculations of the Influence matrices. 

The INFMAN module can therefore handle generally shaped wakes and yet be 
cognizant of the fact that wakes are attached to and dependeun on their genet' 
atlng surfaces. XMFMAN uses the wake corner points generated In POTWAK to 
compute the Influence of these finite vortex segments through the Blot-Savart 
law. This Is then added onto the velocity field computed In VVIM resulting 
from the lifting surfaces to give the resulting perturbation velocity v_. 

Table 9 shows the INPUT deck for the wing-alone case of tables 1, 5, 
and 7. Some of the command formats used are described below: 

1, VNADD — this causes the Influences of one or more wakes to be added to 
the normal Influence matrices (l.e., velocity-dot-normal) that do 
not Include the effects of the shed wakes. 

At this point, the program expects the following data to be entered by the 
namelist $DATA. 

NCOMP — number of lifting surface networks. Default Is 1. 

NIWKS — number of Independent Wake networks. Default Is 1, 

IDGEOM — array of Identification numbers of the geometry files for each 
lifting surface. 

IDWKGM “ corresponding array for the wake file. 

IDIN — array containing the Identification numbers for the normal Influ- 
ence matrices for the lifting surfaces Including the effects of the 
shed wakes. 

XDOUT — array containing the identification numbers for the revised 
matrices Including Che effects of the shed wakes. 

COROPT — integer defining the specific type of core model. Cores ate 
used to alleviate unrealistically high velocities near the finite 
strength line vortices that are used to represent each wake. 

CORPOT ■ 1 creates uniform core sizes. 

COROPT ■ 2 creates core sizes Chat generally depend on Che geometry of 
Che wake. In particular, this option usually results in core sizes 
that are equal to a certain fraction of the "distance*’ between a 
segment and the closest other segment that Is not on the same wake 
line. 

PARAMC ■- Array containing parameters defining the core size. The only 
one of Interest to the user la the first element PARAMC (1), giving 
the core radius corresponding to the COROPT option. 
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2. VEADD this connand causes the Influence of the shed wakes to be added 

to the velocity Influence matrices that do not Include shed wakes. 

3. VIAOD, V2ADD - similar status to VEADD. 

The VNADD command is addressed primarily towards enforcing the boundary 
conditions of equation (22) by forming the \f.N terms. Commands VIADD, V2ADD 
enable the resultant velocity field (v 4* U) to be computed in POTFOR. 

From figure 3, the length of the wing leading edge is 

/3.487A^ + 1 ■ 3.6279. With 16 chordwise panels and for a constant core size 
model COROPT - 1 the panel spacing along the leading edge is given by 
3.6279/16 • 0.2267. Card 5 in table 9 shows the core radius PARAMC(l) to be 
taken as half this value. As before, the geometry and wake files have 
ID « 30701. The ID for the v^.N matrix in card 4 is IDIN ■> 3201, thus over~ 
writing the ID for the corresponding matrix without shed wakes, cards 6, 8, 10 
of table 7. The velocity matrix jv has an Identification IDOUT 3202. 

Table 10 shows the corresponding INFMAN deck for the canard-wing case of 
tables 2, 3, and 8. This consists of two blocks each of which has essentially 
the same !^;itructure as table 9 for the slngle-lif ting-surface case. This uses 
the two geometry files with IDGEOM ■ 8701, 8702. For the canard the wake 
file has IDWKGM ■ 8701. The v.If partitions for the canard corresponding 
to A(.y of equation (23) has IDIN ■ 8801, 3802, thus overwriting the 

corresponding IDIN in table 8 (cards 17, 19, 21, 31, 33, 35). The v matrices 
are stored with new IDOUT ■ 9901, 9902. Similar remarks apply to card 19 of 
table 10 for the wing. The normal Influence matrix of the wing is stored on a 
disk file with IDIN - 8804, 8803, with the velocity influence matrix having 
IDOUT * 9904, 9903. Cards 7 and 20 show that the canard and wing use a uni- 
form core model. 

From figure 5, the length of canard leading edge is 21.685 with 10 span- 
wise panels so that the spacing along the canard leading edge is 2.1685. 

Card 7 shows that the core radius PARAMC(l) has been taken to be half this 
value. Card 20 gives the corresponding value for the core radius associated 
with the wing. It should be remembered that separation from the wing leading 
edge takes place outboard of the canard. 

The INFMAN module consumes the major portion of the CPU time. If the 
program exits on time, restart operations from INFMAN with the commands VNADD, 
VEADD, and VIADD replaced by VI ADD, V2ADD, and STOP, respectively. The only 
cards needing change are card 5 of table 9 and cards 7 and 20 of table 10. 


BCDN Module 

VVIM in conjunction with INFMAN generates the normal influence matrix 
v.£ of equation (22). The BCDN module computes the U.N matrijt, from which 
the singularity strengths can be then calculated by solving the set of linear 
equations (22). These boundary conditions are enforced at one point of each 
panel of the lifting surfaces. 
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Input to BCDN consists of POTFAN geometry flies and user Input. The user 
Input consists of connands and data needed to effect the coionands. Included 
In this data are geometry-file Identification numbers, angles of attack and 
sideslip, and rotation rate vectors. 

Output consists of printed output and boundary-condition files. The 
boundary-condition files are used by PSOLVE (the equation solving program) 
because they form the right sides of the systems of equations that generally 
occur when analyzing aircraft by the panel method. 

Table 11 shows the Input deck for the wing-alone case. Some of the BCDM 
conniands Involved are listed below: 

GREAD — reads In the required geometry Information. For the wing-alone 
case this ID Is 30701. 

BCREAD — reads In boundary condition Information such as ALPHA^angle of 
attack and NSETS«number of sets of boundary conditions to be 
compt;ted. 

CBCV — computes the right-hand side matrix £.N. 

SBCV — stores this in a file; the corresponding ID In table 11 Is 3201. 

Tables 12 and 13 show the BCDN Input structure for the canard-wing case. 
Module BCDN must be called twice, once to set up the canard boundary condi- 
tions using canard geometry file 8701 and store this Information on a file 
with ID ■ 8801. This Is followed by a second call of BCDN for the wing with 
ID X 8702, 8802. For the wing-alone case (table 11) only the angle of attack, 
ALPHA In card 4 needs changing. For the canard-wing case of tables 12 and 13, 
no changes need be made. The angles of attack appearing In cards 4 are sub- 
sequently overwritten In POTFOR, 

The linear equations (22) are solved In PSOLVE for the singularity 
strengths. 


PSOLVE Module 

PSOLVE uses ait LU decomposition method with partitioning developed by the 
Boeing Company. Input to PSOLVE consists of Influence matrix files created by 
WIM and INFMAN, right-hand side vectors created by BCDN, and ut9er Input. The 
user Input consists of commands and data needed to effect the coiomands. 

The output from PSOLVE consists of printed output, solution files, and 
Inverted or partially Inverted Influence matrices. The solution files are 
used by POTFOR to calculate the forces, moments, pressures, etc. 

Tables 14 and 15 display the PSOLVE Input decks for the wing-alone case 
and the canard-wing case, respectively. These decks need not be changed during 
runs. The user Input Is entered by the namelist $DATA and consist of: 
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NCOMF number of lifting surfaces. Default is 1. 

IDIN *- influence matrix ID. 

IDBC — boundary condition ID. 

IDS -* solution file ID created as a result of PSOLVE. 

For the wing-alone case IDIN • 3202 « corresponding to IDOUT ■ 3202, 
card 4 of table 9. Again, IDBC ■ 3201, corresponds to card 7 of table 11, 

The solution ID is now carried over to 3201, overwriting previous Information. 

Table 15 is essentially a generalization of table 14. Elements of the 
IDIN array of table 15 correspond to IDOUT in cards 6, 19 of table 10. 

IDBC * 8801, 8802 is transmitted by cards 7 of tables 12 and 13. The singu- 
larity strengths r^, Xw> corresponding to the canard and wing, respectively, 
are stored in files with IDS ■ 8803, 8804. 


POTFOR Module 

Finally, POTFOR calculates the forces, moments, pressures, etc., using 
the Kutta-Jonkowskl equation 


4jP *■ p Pv X (24) 

Input to POTFOR consists of POTFAN velocity- influence matrix files, solu- 
tion files, geometry files, and user input. The user input consists of com- 
mands and data needed to effect the commands. The seven major POTFOR commands 
are the following (tables 16 and 17). 

1 . READ reads in geometry and singularity data for a network. The follow- 

ing data is expected from the $DATA namelist. 

IDG “ ID of geometry file. 

Foe table 16, this is 30701, while for the canard-wing case of table 17 
IDG - 87Q1 for the canard file and IDG ■ 8702 for the wing geometry file 
(cards 4, 5, 31 and 32). 

2. ATTITUDE — this causes attitude data to be read in. The following data 

is entered via $DATA namelist. 

IDSF ■ ID of solution file. This is the output from PSOLVE 
module . 

For the wing-alone case IDSF ■ 3201, card 3 of table 14. For the canard-wing 
case IDSF ■ 8803, 8804, from card 2 of table 15. 

3. COMBOS — this command determines a transformation matrix from which the 

"combinations" may be determined from the "cases." By definition, a 
case corresponds to a boundary condition vector calculated by the 
boundary condition program and for which a solution vector was 
determined by the equation-solving program. A combination is a 
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linear combination of cases. For example, the solutions for two 
different angles of attack may be combined to yield solutions for 
many different angles of attack. 

This command must be preceded by a READ command. All groups are dropped 
except those containing the geometry information. Subsequent comnands may 
thus overlay whatever force information has been computed at this point. 

The COMBOS command Is reqv‘ red after every ATTITUDE command provided that 
the desired combinations are different from the cases. 

The program expects the following Input by the $DATA namelist: 

NCOMB " number of combinations. 

ALPHA - angle of attack. 

A. NETLOADS — this calculates the net forces and moments using the Kutta> 

Jonkowskl equation (2A) . The following data Is Input from the $DATA 
namelist: 

NCOMF "< total number of lifting surfaces. Default Is 1. 

IDVE ■ identification numbers for the boundary condition veloc- 
ity influ<*,nce matrices. 

IDS ID for the solution vector of the lifting surfaces that 
have a velocity Influence on the current surface. 

PRINTV >■ logical variable whose truth causes the net velocities 
at the boundary condition points to be printed. Default 
Is .FALSE. 

In table 16, IDVE corresponds to IDOUT *• 3202 of the INFMAN module 
(table 9), while IDS corresponds to IDS ■ 3201 of PSOLVE (table 14). A 
similar generalization holds for table 17, In comparison with tables 10 and 15. 

5. PRESSURES — computes pressure distributions on the lifting surfaces. 

This requires the following data: 

POPT ■ 0 computes pressure In the usual way (A£.N/area) , the 
default value. 

6. SPANLOADS — calculates span loads and requires the following Input In 

$DATA namelist: 

SOPT - Integer variable whose value determines what types of 
spanloads will be calculated. Default Is 0. This causes 
the spanloads to be normalized by the local chord and 
resolves in the wind-centered coordinates. This produces 
section lift and Induced drag. 
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SOFT ■ 1 spanloads normallzad by avaraga chord CAVG. 

CAVG - FLT(3) of POTGEM 

7. NETK — similar in purpose and organization to the NETLOADS conmiand, but 
can also Include edge suction forces for thin lifting surfaces. 

From the user standpoint however, table 16 for the wlng^alone case needs 
no changes, while table 17 for the canard-wing case requires only one change 
for the angle of attack at card 9. 


POTFOR Output 

The POTFOR output for the wing-alone case is displayed in table 18. The 
forces, moments, and pressures are calculated from equation (24). The POTFAR 
program uses two versions of this formula: an approximate or modified version 

that uses the mean velocity over each panel. This is computationally faster 
and better able to model curved surfaces. The exact version uses the correct 
velocity V on the panel boundary and gives good agreement with experiment 
for the wing-alone case. The approximate version consistently gave worse 
results. However, for the canard-wing case, sometimes the approximate ver- 
sion performed better. Xn table 18, the first eight outputs correspond to the 
approximate version and the last three to the exact version. 

1. VELOCITY AT CONTROL POINT LOCATIONS; This contains the velocity compo- 

nents u, V, and w at each of the panel control points at which the 
boundary conditions 22 are enforced. 13 ■ 1 indicates the wing has 
no thickness or camber. 

The index 12 runs from 1 to NBPV - 8, which is the number of spanwise panels 
and II runs from 1 to NBPS • 16, the number of chordwlse panels. The ray 
nearest the leading edge corresponds to 12 ■ 1. The u,v components on this 
ray are fairly uniform near the apex and trailing edge, indicating nearly 
conical flow. However, within the Interior of the wing, the flow field is 
less conical. The w-components throughout are very small, as Indeed they 
must be, because the flow must be tangential to the wing as Indicated by equa- 
tion (22) . 

2. NORMF, NORMM, ADD: These are defined by: 

NORMF > normalization factor to convert dimensional forces to 
coefficients; 

NORMM(l,2,3) * normalization factors to convert dimensional 
moments to nondimens ional coefficients; 

ADD - logical variable whose truth causes accumulation of 
forces and moments. 
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NORMF, NORMM(l,2,3) arc defined by: 


NORMF - 1/FLT(1); 

NORMM(l) - NORMF/FLT(2); 

NORMM(2) - NORMF/FLT(3); 

NORMM(3) - NORMF/FLT(4) 

where FLT Is accessed through the geometry files. Table 1 shows that 

FLT ■ (1.7437f 1.0» 3.4874, 1.0), so that NORMI' ■ 0.573 as Indicated In the 

output display. 

3. COMBINATION;! lifting surface. (FX,FY,FZ) and (M0MX,M0MY,M0MZ) are the 

net force and moment coefficients on the network. (FXW,F^,FZW) are 
the force coefficients In the wing centered coordinate system with 
FXW - drag, FYW - side force and FZW - lift. (M0MXW,M0MYW,M0MZW) 
are corresponding the moments coefficients with MOMXW " rolling 
mo: t, MOMYW ■ pitching moment, and MOMZW ■ yaw moment. 

4. COMPUTATION OF PRESSURES; The default value of SOPT in the PRESSURE 

command Is POPT * 0. In figure 3, S,V corresponds to the x,y 
coordinates. Table 18 then ohows the spanwlse variation of pressure 
coefficients. 

5. SPANWISE LOAD DISTRIBUTION; Reference chord - CAVG - FLT(3) of POTGEM. 

Components of the force on each spanwlse row of panels are FXBW, 
FYBW,FZBW. These forces are made nondlmenslonal by the product of 
the dynamic pressure and CAVG. These are expressed In the wind cen- 
tered coordinate system. For a wing with ray paneling the longi- 
tudinal load coefficient is given by 


C 

n 


I. 


-s 


b/2 


y^XBW)2 + (FZBW)^ 


(25) 


where b Is the wing span and s is the local semispan. 

6. VELOCITY AT NI FORCE SENSING LOCATIONS 

7. NET FORCES CAUSD) BY NI VORTEX SEGMENTS 

8. VELOCITY AT N2 FORCE SENSING LOCATIONS; This velocity Information Is not 

directly relevant to the user. It contains Intermediate calcula- 
tions for assembling the velocity contribution In the NI, N2 direc- 
tions. The correct use of equation (24) requires the velocities be 
calculated along the panel boundaries rather than using a mean veloc- 
ity at the control points. POTFAN arranges the corner points in a 
rectangular grid spanned by directions NI, N2 corresponding the 
number of corner points along the leading edge, and number of points 
characterizing each shed wake. 
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The compoelte effect ie described In the following heedingSy which give 
the output of the second method of computing the forces end pressures. 

9. NET FORCES CAUSED BY N1 AND N2 VORTEX SEGMENTS: These give the force end 

moment coefficients for the exect use of 24. 

10. COMPUTATION OF PRESSURE 

11. SPANUISE LOAD DISTRIBUTIONS 

The leyout of the output for the canerd-wing cese is virtuelly the seme, 
and comes in two blocks corresponding to the canard and wing, respectively. 


The Advanced Panel Code 

The user's manual for operating the pilot version of the APC is given in 
reference 10. The user has to prepare a subroutine INPUT written according to 
the specifications given in reference 10. This must contain in particular, a 
mechanism for calculating and accessing the wake corner points for the leading- 
edge vortex trajectories. This interfacing is done through POTFAN geometry 
and wake files. 

Table 19 contains a listing of the subroutine INPUT for the wing-alone 
case already encountered in POTFAN. The user input is entered through name- 
lists DATAO, . . .,DATA4. The variables involved are defined through a series 
of comment cards. The program distinguishes lifting surface networks and wake 
networks through an integer NTD. which is 18 for a wake and 12 for a lifting 
surface. 

There are two wake networks: a wake network from the leading edge and a 

semi-infinite sheet from the trailing edge. The mst convenient way of gen- 
erating the latter is to create a second fictitious semi-infinite rectangular 
wing attached to the trailing edge of the actual wing. This is done through a 
POTGEM file shown in table 20, which is very similar to table 1. 

A rectangular wing of length 9500 units and semispan 1 is used. This 
arrangement has 15 panels spanwise and one chordwise panel. The ROSS command 
is used in cards 30, 31 to position the wake at the wing trailing edge. The 
POTGEM file for the actual wing is Identical to that in table 1. While POTFAN 
was not sensitive to the direction of the final seml-lnf lalte vortex segment 
from the leading edge, the AFC was sensitive to this direction. It was found 
desirable to orient this line at an angle a/2 to the plane of the wing. So 
that FLT(5) on card 33 must be replaced by cos a/2, 0, sin a/2. It should be 
noted that NBPS for the traillng-edge wake, card 6 of table 20 must have the 
same value of NBFV for the wing, so that the grid point along a common bound- 
ary are identical. 

The data for the subroutine INPUT is given in table 21; almost all of 
this data block are APC variables which are described in ref e;i ence 10. The 
data organization is summarized below: 
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First rssd ststmnsnt mt 11ns 71* 

Second read ststement at line 218. Read in $DATA2 namelist. 

Third read statement at line 221. Read in $DATA3, within the DOLOOPIOOO. 

Fourth read statement at line 2S7 $DATA1 

Fifth read statement at line 259 $DATA4 

Sixth read statement at line 362 $DATA0 

Seventh read statement at line 381 — $DATA1 

Eighth read statement at line 383 ~ $DATA4 

The leading-edge wake corner points are accessed between lines 360 and 
379. This is the information contained in table 6. The subroutines OPEHR 
and RELFIt are part of the POTGEM module and are described in reference 6. 

The APC uses a main overlay program called FEE, wh^ch must have its argument 
list expanded up to tape 20 to be compatible with the POTFAN disk file system 
so that 

OVERLAY (FEE. 0,0) 

PROGRAM fee (INPUT, OUTPUT, TAPEl, 

TAPE 4, TAPE 5 - INPUT, TAPE 6 - OUTPUT, 

TAPE 7 TAPE 20) 


• 

STOP 

END 

This sequence of cards is placed at the end of the subroutine INPUT. The 
job control cards are shown in table 22. The program organization is shown in 
figure 8. 

A major drawback of the technique is the current upper limit of a 1000 
panels for all the networks. To get realistic results for the wing-alone case 
requires about 990 panels. It is felt that currently comparable results are 
not possible for the canard-wing case with this restriction. 
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OF POOR QUALITY 


TABLE 2.- FOTFAN GEOMETRY INPUT FOR CANARD (TRAPEZOIDAL PLANFORM, STANDARD 

PANELS) 
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S.t .gv»^ iwf nl«Trn''NS ii» Tm, ••«(<(> 

51 • I -I 

S,» • 1 .1 

Vi • (1.0 

5V • 1 1 

V5 • •» •! 

'jfc » (i,g ,J«?Ur 

5T *0 

5(* •MMIbn 

50 • Ivt(1(ij^i, FL ri5t*. 1,^5, 11,51,31, is, a, ri r(io)i%^,|o teNii 

i>() •SU’xt 
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hi atiVIHT 
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TABLE 3.- POTFAN GEOMETRY INPUT FOR WING (TRAPEZOIDAL PLAKFORM, STANDARD 

PANELS) 


I 

« 

I 

« 

1 

It 

4 

10 

u 

IJ 

l» 

n 

10 
ir 

It 

!• 

20 

It 

11 
ii 
l« 
l(i 
I* 
IT 
Id 
14 

in 

it 

il 

Si 

id 

s^ 

in 

iT 

id 

14 

40 
“I 

41 
4 i 
4« 
4 b 
«l| 
4T 


• I N »• U T 

•flfti * ' * * ** 

*?N o-Ftiod N|Nfi r iNP.Lftt M4itr.(MlNa sinttjTiiv nut, 
•c*»»r 

•OSlCdlNti 

• t04M NHnilT,NttrVi4 I£N|> 

• SHIt 

• iiNC*yi |C4|,nidTto itNf. 

4 fiNCdVi |C4|,Cn»iil »Chn 

• SHlI 

• 4INC441 |(■■H,c(1|•r«4,V4ll|••40 |luU 
•9«n 

« H4C4V) Ic45,c0*>ti0 IJNh 

• lINCttVt tr»T,C0dTl4,V4KI>l «f»)n 

• SmuntflNt NINO 19 |I4 fl<r lift OiSnI) 

• ilNC4V| ?f*(i,cndliO ifNP 

•)ti> (itinNt »Mt »nnf) 

• tOill (lldT9Vi«,viH|SV4i,d| IfNO 

•sn totriHt thp locifinw od rxir rtiiTHitittj cinmp vonnt 

• tOlTl 41«ji9V»|4,itO 4l«P 
•vL (iitnNr t«e (.riolNc lOtti 

• IDITI tn»T9Vii|,4flMSVi|,Vlt»(),V4,,n,|1,40. 

. vi8l3V44T.«5,«b.ftM ifNi) 

• VU (DCrtNf TNC THIII 140 CItbt ) 

• tliiti vm?9VitT,is,4(,M limn 

• ViMC (OtrtNt TWr 981N4t»r "INrl. r,|9T»tMI)M0mt 

• iPiTi iii(*t» 4 limn 

•BLUC (ittrtNt rnr c'tOBnttist dimiL niiTOlHiM inNi 

• »o*u iiit*T4i irwo 

• ivic (injwsT tMf •imets so tmt ihcv dircM c*Ni«ut 

• lOltl )l*nJi,|li«i5«4T4,Nlt)JlO »r40 

»iviic (io,|U!ir iHf diMfiS 90 ?Mir TMf,r tircM cini»o) 

• lOiti itiin 

• (MHInI Ol)t 5>44iL9) 

• 4(id,4|4fL Urn rHt cimiMu 

• mtNi 

• lOMi 4sii» irmo 

*B099..M0|| ClNlBO TO LiTl tlY 1,»| ‘(mrtS 

• |01U 4H|tO,fl,<l!|MtniO,0»»l,(>t fU* 0 Il4rj 

•09M. toesjNi iH( oouoirt stmoijiltiTT rufisi 

*1-1 

• I -t 

• i lmrr»tmt si4oui.ioIt» 

• ,M! -1 

• 4 i.r‘0140 roul stmsui «nTT 

• I ♦! 


44 

tin 

SI 

Si 

51 

Sil 

ss 

94 

9T 

54 

54 

no 

M 

III 

41 

44 


45 

4b 

il 

44 

44 


• -I •! 

• 4 tmiuso tiifit 5rtr.MUi«rTT 

*•1-1 

• 1 

• IT ifi’ simoiiuouT 

*n 

•t‘v» (iteKimi ihf niscrnoss nr inr sino «*4rst 

• I -I 

• 1 *1 

• .41404 0,n ,14101 

• II 

• -I -I 

• .oiono n,n ,la?ul 

•0 

•rimriiH 

• K'lll ISKifttni, riTtSIb. 1,15.4, 51.1), 15. a, rlt(|nt«S4.|U lisu 

• STI’44 

I iu«i7np »rst' 

• pi»|*ll 

« null »f„n 

• STOP 
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TABLE A.- PARTIAL LISTING OF CANARD POTGEM FILE FOR INPUT OF TABLE 2 


;imi 

tN O-firtS II HiftC 


.oir6Mt*trs 

»SI*II 

rS«t) 

♦ SRI) 

^SRIt 

♦ Silll 

i*»n(oeR»Ne ming i*( ihc i«o 
| sL (oifiwe t«r 

^svt tORriNK ThI unc*f|flM Or th* tiri 

•Vi. rncriNt rne irtniMr. eociri 

;vu (oeMNi fMe immiins enon 

fVtHC (Gtiti^jt iMf nis?*i«yTiGN) 

♦ StGC iPtRlNt I«r CHn*nR|U l**NeU OI<|T»tH'inON| 

1 HO“»R»Rtl. U? lM|t C*NtRp 

*H0»S—»0Ve e»»*»"n to Irfl sr »,*l iimra 
*osrt inmNi iHf omutri stNcMi.»fiTy Fi*r.sj 

'.uvt« inf F INC fHt nlRf.CMONS OF IHk RHfU a*FMI 

iFlNfSH 

.jfORF 

t.iNlT I 0F«NkO Fn» RRITIMG FU,e OF Tv»»e I 4NO IGa efPI 
5«e*MIN l|t<f ■ On/j|/M lJ.JI.0f 

T*Ff 1 H*J hUN fROFII.eO« STlIRkD »» rILr FF« f 0 IGHf « *"0 «l*<OUhOf 
r»r*T|rtN (IF GFlfMkTRV Fitr 


tlTLE I IN n»H(Ik f*N4H0 H iNClfO HfUr 


fUOG) » fFlFtFFFFFFfFFFFTflF 

}INf) • OtlklOktOOOIOOO 
(FUl ■ Jlk.lOAAll 25,400000 >1,310000 25,400000 1,0000000 0, 

.12127*700.01 50 . 1*0000 0, 

♦P*|N| 

pHINfnni OF GtONM** Mif 0414 


rliLf • IN n.i505 CiNfon II 4 NCU 0 *4><e 


r4C*T|0N Hie • 11.21.07 

ii^op'*) » uioiiiim 

flO) * tifni 

aOU) ■ friFTFFFFFFFFFr 
flNIl 1 0 h * 10 5 0 

fFLTi I 5l«,inoonno ?5,4o«00on 
121,278ft45? .0321274 


0 0 0 
21.1100000 
5*. I’lOoooo 


I 0 II 0 

2S,«oooonu l.nonoonii u.ooooooii 

n.ooooonu o.ouoouno 


321^2fNfu 


O.OIIflPOOU 


’oFNkl, CnUNkN P'llNfS ; 

,.__j — rrTTTT*” TM.J) 

1 1 io.'»2ii40oo fl.onnooog 

2 1 12, 512*07* 1.1t!2l»S 


f(|..l) 

o.ooonnoo 

o.niinuouo 


5(Ii Jl 

l,*l oduoii 

5, 1212105 


V<l, J» 
10,4205000 
iJ.SIPlRf* 


1 

1 

1 U. 

2.0224 J4(i 

o,flO«rtOOi' 

B. 4124100 

1«. ||)4*4S| 

u 

1 

iS.dOofO^l 

l.'MlbSOS 

n.nuoniion 

7,7*3i»5«5 

15. *4*7*27 

5 

1 

J7 .2N4««I/? 

5.244*7*0 

O.OOOOOOO 

4,gS4|)7*t1 

17,2**11*112 

h 

1 

JO.ONOO*!* 

B , ^HortoTo 

(1.111)0(10011 

lO.lNBORf* 

1*. 0*04*7* 

7 

1 

J0.47 1(,45« 

7, *071171 

O.OO'IOOOO 

1 1 ,077 11 7 1 

40,4710*54 

5 

1 

42.00*1*24 

4,1 7*5l4b 

0,(1000000 

12,4**51*0 

42,0*51*24 

4 

1 

>il.057?*n5 

10,4*47501 

0,0000000 

14.24475*1 

4l,**72SV5 

10 



1 1 ,*ftt»47S«> 

o.ooonooo 

15. *11)475* 

45,24417*0 

1 1 

1 

OB.O'I 1 i| fS* 

11,1121451 

o.ooooeuo 

1 *. 4221451 

4*. 041*75* 

1 

2 

17, ■•42*00(1 

11 , 6 oim ffii'o 

O.OOOOOrtO 

1.410000b 

14, 04250(10 

2 

2 

}S , to 1 n loo 

1 .11121 45 

0,0000000 

5.1212145 

15,76101** 

1 

2 

17.074*712 

2.O224J40 

O.OOOOOOO 

*,4l24l*V 

17.11745 7 12 

u 

? 

1*, jn»,no« 

1,411*5*5 

0 . onooono 

7,701*5*5 

10,)4M|U4N 

5 

2 

1 ’’, 7lB„Uft* 

5,244*7*0 

0 , 000000(1 

4,iI54«7*u 

14 , 7|*U<I*1 

«r 

' 

>1 > 1* . •* JO 

. ««Oll4Tf. 

0 oonooiio 

10.1**0476 

41.0151*24 
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TABLE 5.- POTWAK INPUT DECK FOR DELTA WING 


$DATA 

NSEGS 

m 

X50, 

ALPHA 

n 

f2E+02, 

MACH 

m 

0.0, 

UINF 

m 

.XE+OX, 

ID 

m 

3070Xf 

PLOT 

m 

T, 

STORE 

m 

T, 

IDWAKE 

m 

3070X, 

DELTAT 

- 

,XE+Q0, 

RX 

at 

.X3429E-f00, 

R2 

m 

.69487E-0X, 

THETA 

» 

,95833E+00f 

GAMO 

m. 

.22767E+00, 

STASHN 

m 

X, 

FACTR 

M 

.5E+00, .XE+OX 

XV AL 

m 

.XE+OXf 

DELTAX 

m 

.3E-0X, 

EPS 

m 

. XE-OX, 

IPLAG 

m 

0, 

ILINE 

m 

2, 

IVLO 

m 

X25, 

lYHI 


925, 

IZLO 

m 

X25, 

IZHI 

** 

925, 

lARGMIN 

M 

X25, 

IARGM'aX 

9B 

925, 

IRKN 

SC 

xo, 

VLB 

M 

.XE-02, .XE-02 

VUB 

M 

.SE+03, .5E+03 


“Ij “If “If 


lE+OX, .lE+Ol, 


0.0, -.5E+03, -1, - 
.1S707963268E+01, . 
“X 


1# “If “If “X, 
5E+03, 
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IPRINT 

» 

1, 

NCON 

m 

1, 

NSIDE 

m 

1, 

INFO 

m 

1, 

ITMAX 

m 

500, 

ICNDIR 

» 

0, 

NSCAL 

m 

0, 

SCAL 

m 

-1. -If 

ITRM 

m 

2, 

CTL 

m 

lE-01 

CTLMIN 

s 

CM 

O 

1 

H 

PHE 

- 

O 

• 

o 

DELPUN 

m 

.lE-08, 

DABFUN 

- 

.lE-08, 

FDCH 

- 

.lE-09, 

FDCHM 

m 

.lE-14, 

THE 

m. 

o 

• 

o 

NFDG 

m 

Of 

LINOBJ 

a 

0, 

CT 

a 

-.5E-01 

CTMIN 

$END 

a 

o 

o 


CONCLUDED. 
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TABLE 6.- WAKE GEOMETRY OUTPUT FILES 


1ST RECORD 

NTITL=40 

CTITLCNTITL))--ALPHANUMERIC TITLING INFORMATION FROM THE GEOMETRY 
FILE CORRESPONDING TO THE COMPONENT TO WHlCH THE WAKE IS 
ATTACHED. 

NRECS=5 

CIFORMCNRECS)) — I, 1, 0, 0, 1 

NID=2 

!DC1)-~ID NUMBER OF SURFACE TO WHICH THE VORTEX IS ATTACHED. 

ID(2)— ID NUMBER OF WAKE FILE 

NL0G=1 

LOGO)— NOT USED. 

NINT=4 

INTCO--NSEGS, THE NUMBER OF CORNER POINTS AND SEGMENTS DESCRIBING 
EACH SHED VORTEX LINE. THE FINAL SEGMENT IS AN INFINITE 
LENGTH SEGMENT DEFINED PARTLY BY DIRECTION VECTORS. THUS 
THE NUMBER OF CORNER POINTS IS REALLY NSEGS AND NOT NSEGS-1. 

INTC2)--NLINES, THE NUMBER OF SHED VORTEX LINES. 

1NTC3)— Nl, THE NUMBER OF Nl DIRECTION CORNER POINTS ON THE 
COMPONENT FROM WHICH THE WAKE IS SHED. 

INTC4)--N2, THE NUMBER OF N2 DIRECTION CORNER POINTS ON THE 
COMPONENT FROM WHICH THE WAKE IS SHED. 

NFLT=3 

FLTCD--NOT USED, BUT RESERVED FOR MACH NUMBER, 

FLTC2)--ANGLE OF ATTACK CORRESPONDING TO THE WAKE GEOMETRY, 

FLTC3)— UINF, THE FREE STREAM VELOCITY. USUALLY JUST 1.0. 
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TABLE 6.- CONTINUED. 


RECORD 2 
JUNSEGS 
J2=NLINES 
J3=3 

NW=NSEGS*NLINES»3 

CXWAKECNSEGS, NUNES)), CYWAKECNSEGS, NUNES)), CZWAKECNSEGS,NL INES))- 
X, Y, AND 2 POINTS DEFINING, IN PART, THE GEOMETRY OF THE SHED 
WAKE. THE WAKE IS COMPOSED OF NUNES LINES AND EACH LINE IS 
REPRESENTED BY NSEGS-l FINITE SEGMENTS AND ONE SEMI INFINITE 
SEGMENT. 

RECORD 3 
J1=NLINES 
J2 = l 
J3 = l 

NW=NLINES 

C I 1 INDX(NLINES))— ARRAY GIVING THE II INDEX OF THE CORNER POINT 
ON THE GENERATING SURFACE FROM WHICH THE WAKE IS SHED. 

RECORD 4 

J1=NLINES 
J2=l 
J3 = l 

NW=NLINES 

CI2INDX(NLINES))— ARRAY GIVING THE 12 INDEX OF THE CORNER POINT 
ON THE GENERATING SURFACE FROM WHICH THE WAKE IS SHED. FOR 
EXAMPLE, THE 4TH VORTEX LINE ORIGINATES AT THE I1=I1INDXC4) 

AND I2=I2INDXC4) CORNER POINT CORRESPONDING TO THE GEOMETRY 
FILE WHOSE IDENTIFICATION NUMBER IS 10(1). 
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TABLE 6.- CONCLUDED. 


RECORD 5 
Ul=NLINES 
J2 = 3 
J3=l 

NW=3*NUNES 

CUVWXCNLINES)), CUVWYCNLINES)), (uvwz(nlines))--arrays that define 
THE DIRECTION OF THE FINAL SEGMENT IN EACH VORTEX LINE. NOTE 
THAT THE FINAL SEGMENT IN EACH LINE IS INFINITE IN LENGTH. 
ALSO, BY DEFINITION, UVWX( I)**2+UVWY( I )**2+UVWZ( I )#*2=l 
WHERE 1 .LE. I ,LE. NLINES. 
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TABL.'? 7.- WIM INPUT DECK FOR DELTA WING 


■a i itM Mfla m u a^u m ti Mt i i u i«» 

’ wniioi 

’ TT»r'r‘»5,»a*/T« ts'.an.i* 


t •!■ »* It T 


2 • IO*Ti UNO 

I — — — — ' ■ 

• • tO*T4 lOIfNOaiOTOI ICNO 

% •44w»«c4UUi'4u iNfiuCMH « m<Kifr »f n.c. noinis 

• • |0*T4 IO*iaoi,N«HNtkal UNO 

r *jim-c4Ui4«i4 4» N| roNcc t.nc*Mo*>i 

• • I0»l* IDiiaOl UNO 

to • »0*t* I0»l20l ICNO 
II 1U«* 


UM4I 4 NO«^ |N«M MLC 
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•oiris voMK ixriitMcf i.* 


nV'<*'*ic MMUttyi }AAM wnini 

Ti"» ■ 11. ?*,(?• 


« I s * <• t 


i •0*1* 

? • fOtU NRiiM«yi| tUNP 

> • 

« * Hi Ij* 1n*i,r RiCToM 

I • I* ■■■■■■■•til 

t 

T • io»f* «cRn 
R • 

f INftUrNCl Pf iHf. C*NR*0 KH iTitlR 

11 • 

12 •RH(«0 

II t toil* |p$fr>m«i|7|U »RNP 
m tcnHi* 

It • tO*M RAChI*.} URO 
Ik •!Nl| 

i; • io«it litiaMflt «rNn 

IR •viil 

I* • lOti* lo««)ini «2Nn 

20 •VIJ2 

21 • IIXI* IlIkkRAI iFNn 

22 • 

2S • tOHRUTr I'truiruct nr ipr e»N»ki) on thi <i|Ns 

2>i • 

25 * 

2 k «nofti) 

27 • «DM« lllRf riRTO? tlNO 

2 R •c'OMk 

2k • IIUI* ItHn 

III •I'J^l 

II • ll)M« ll>kkR02 (FRO 
J2 *VIM 

I) • titiU in>iiR02 |FMI> 

U 'ylSF 

IS • «li«U ltl■*>IU2 |R»in 

Ik • 

I; • ci»«i»uTr inriiiysfi nr iKf k|xri u»i iMiiF 

IK • »•••••, •«n»mi».»»*kikk«, Ik, ■■•»kkki« 

19 • 

tin kiiRCui 

'll • <u>t* iiiKfNnikriiP 17 IIP 
112 

»1 • VtIM 

<lk klRM 

ns • tn«i4 liiikHni tn'iQ 
■Ik •VISI 

u? • tntit iiikiiHni kfNR 
1 C •vlS2 


■Ik • t)|k|* UKkNUI ^fNfl 

5(1 • 

51 • cimifr I'lriiirnri or iiif unr ji »hr r»k|ko 

52 * mkkkvkk, 

SJ . • * 

5u •kiltkU 

5s • i'ui* lokfrik/oi Hip 

Se •(•U'p 

57 • iPktk tiMn 

5n • I 111 

51 • 7P«t» IPikCna trip 

kO »yl!i| 

M • It'll* IlikMail tf’iP 
k2 *VI52 

II • ii'iii iPikMOk kr’ip 
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TABLE 9.- INEWAN INPUT DECK FOR DELTA WING 

lai i m lu u u i umm 4< n »M ■ m uta, • 

ovmamic mcmonv • «|••^wa•o■ 

t* - ■ • 

Tp!rTli/i»7fi uToTiii’ ™ 


— I 


.4- 

2 
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• |0|Hl|2t|»|OoU!iSI02. 
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k •VCtOD 
I ■ UNO 

• tVUOD 
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TABLE 10,- INimK INPUT DECK FOR CANARD-WING CASE 




nv****ic I Mjnno tto»nii 




r H ^ II ! 


I 

•> 

I 

!l 

« 

k 

T 

» 

k 

10 

11 

12 
1} 
l« 

Ik 

IT 
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lO 
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21 

22 

ti 

iu 

2« 

2k 
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TABLE 11.- BCDN INPUT DECK FOR WING-ALONE CASE 


Mouum. lauMfliiay c owti fn n viMioti 


»y«ii*ic Mctosy ■ !«•«• «o*oi 




aaaaaaaaaaa aaaaaaaaa- aaaaaaaaa •' 
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J 1 - 

) •ICMAO 
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7 • )|tl 
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TABLE 12,- BCDN INPUT DECK FOR CANARD-WING CASE (CANARD) 


nniwi>»<<* epNpjTfoN •nnoi»*H, vrutfuN i,; 


(>*N»M|e • lAfloA 






I N » ti I 
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? 

• nnii 



1 




11 


n|*N*al«,Ntf.Tli< IFNI) 


H 
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1 
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in 
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II 
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1? 

•»Tim 




lINft It fS nOh twr l*<»IIT MLt 


original Pf|SE b 

OF POOR QUALITf 


■ M . 


TABLE 13.- BCDN INPUT DECK FOR CANARD-WING CASE (WING) 


•ofpAN cn"PirinN R#no»*"« vr*»*fO»* i#i 

nfXAmC '§€»‘liM* ■ lo<i5« MRsps 


tl*«l « M. 111 . ft* 



• ( V M II r 


1 

*r.»l*u In r.enNfTiiY uf iwr i.|j6i 


i 

• ftTOi 


\ 



tf 

* «»CN6*0 (♦•Nft 


s 

•CftCV 


* 

• 9 *Cv 


T 



A 

•NC»f*0 


4 

♦ »«c*E*o »rNft 


ID 

• C*Cv 


II 

• 5*Cv 


i; 

• STOP 
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TABLE 14.- PSOLVE INPUT DECK FOR WING-ALONE CASE 



(LU DtCOHfOllliaN l.|TH0n UIINS PtHf MIONtHd AND »*Nt)nM )/0> 
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TirT Of' »u»' u-^ ' '*’“ - ■» 
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(a'' Coordinates for cross section ' I 

(ref. 1). I 

Figure 1.- Mangler-Sraith Vortex K. Mup Analysis. ' 1 

I 

'i 
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(b) Approximation sheet 
(ref. 1). 

Figure 1.- Continued. 
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(c) Transform of figure .l<b) . 
Figure 1.- Concluded. 
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Flgur* 2.- POTFAN system. 
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Figure 3.- Singularity panels with flow separation on all sides. 
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Figur# 4.- CloM-coupl«d canard wing-«odal of rafaranca 8. 
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Figure 5.- Planform boundary Curves In POTGEM. 
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Figure 6.- Vortex rollup for delta wing at 20*^ angle of attack. 
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TOP VIEW 



SIDE VIEW 



FRONT VIEW 



Figure 7.- Different views of trajectories. 
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Figure 8.- Program structure for APC with POTGEM.POTWAK, 
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